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Magneto-optical trapping forces for molecules are far weaker than for alkali atoms because the 
photon scattering rate is reduced when there are multiple ground states, and because of optical 
pumping into dark states. The force is further reduced when the upper state has a much smaller 
Zeeman splitting than the lower state. We use a rate model to estimate the strength of the trapping 
and damping forces in a magneto-optical trap (MOT) of CaF molecules, using either the A^ni /2 
- transition or the B^E+ - X^E^ transition. We identify a new mechanism of magneto¬ 

optical trapping that arises when, in each beam of the MOT, two laser components with opposite 
polarizations and different detunings address the same transition. This mechanism produces a strong 
trapping force even when the upper state has little or no Zeeman splitting. It is the main mechanism 
responsible for the trapping force when the A^ni /2 - X^E"*" transition is used. 

PACS numbers: 37.10.Mn, 37.10.Pq, 37.10.Vz 


I. INTRODUCTION 

There is currently great interest in cooling molecules 
to very low temperatures, motivated by a diverse range 
of applications [T]. Laser cooling, which has been used 
to cool atoms to ultracold temperatures for decades, is 
difficult to apply to molecules because it is necessary to 
address multiple vibrational branches, each requiring a 
separate laser. Despite this difficulty, laser cooling has 
been demonstrated for the diatomic radicals SrF [a [3], 
YO d] and CaF [3, and most recently a magneto-optical 
trap (MOT) of SrF molecules was demonstrated [51 [7] ■ 
In all these cases, the laser cooling transition was the 
A^IVi /2 - transition. The mechanism of magneto¬ 

optical trapping for these diatomic molecules has been 
elucidated in [5| where rate equations are used to model 
the interaction of the multi-level molecules with multiple 
frequencies of laser light. 

For a two-level atom of mass M, wavelength A and 
decay rate F, the maximum achievable acceleration is 
hV/{2MX). Using the parameters of the A^ni /2 - X^E+ 
transition of CaF, this is 290,000 ms“^. For the real 
molecule the acceleration is greatly reduced. This is 
partly because the molecule has multiple levels in the 
ground state that all need to be driven with separate 
laser frequencies, and this reduces the maximum achiev¬ 
able scattering rate relative to that of a two-level system 
[9|. This is compounded by optical pumping into states 
that are dark to the polarization of the laser beam that 
pushes displaced molecules back to the centre. Finally, 
even in the absence of these other difficulties, the trap¬ 
ping forces are very weak when the upper state g-factor, 
gu, is much smaller than the lower state g-factor, gi [5]. 
In the Hund’s case (a) limit a state has a negligibly 
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small g-factor and so is unsuitable for making a MOT. 
The observed MOT of SrF was attributed to the mixing 
of the A^ni /2 state with the state, resulting in a 

value of gu that is large enough for weak trapping |8]. 
Magneto-optical trapping of CaF is being pursued, but 
given the very small value of in the A^IVi /2 state [10] it 
is unclear whether the A^ni /2 - X^S+ transition is suit¬ 
able for the MOT, or whether it is better to use the B^E+ 
- X^E"*" transition. In this paper we explore this ques¬ 
tion and we identify a new mechanism of magneto-optical 
trapping that can give strong trapping forces irrespec¬ 
tive of the upper-state ^-factor, and that also side-steps 
the problem of optical pumping into dark states. We fo¬ 
cus on CaF, though our conclusions also apply to other 
molecules of interest. 


II. MODEL 

We model magneto-optical trapping of CaF using a 
rate equation approach identical to that of [8], except 
that here we use the full non-linear Zeeman shifts dis¬ 
cussed below, rather than assuming linear shifts. We also 
use the same notation as in [5|. The angular momenta 
of the lower and upper states are written as F) and Fu 
respectively, their projections onto the x-axis are Mi and 
Mu, and in the linear Zeeman shift approximation the 
(/-factors are gi and gu- Along each coordinate axis of 
the MOT, molecules interact with counter-propagating 
beams, one pushing towards the centre (the restoring 
beam) and the other away from the centre (the anti¬ 
restoring beam). In a local coordinate system, which is 
always aligned to the direction of the magnetic held at the 
position of the molecule, the polarization of the restoring 
beam is (t± according to whether it excites AMp = ±1 
transitions. Throughout the discussion we take the held 
gradient in the z-direction, which is along the axis of the 
magnetic held coils, to be 20G/cm. The intensity distri- 
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bution of the laser beams is Gaussian with a 1/e^ radius 
of 12 mm. 

A useful way to summarize the trapping force of the 
MOT is to calculate the acceleration of a stationary 
molecule versus its displacement along the z-axis, and 
the acceleration of a molecule at the centre of the MOT 
versus its speed in the z-direction. We call these ac¬ 
celeration curves. For small values of the displace¬ 
ment, z, and speed, Uz, we can write the acceleration 
as 02 = —Lu'^z — (3vz, where uizH^ir) is the trap frequency 
and /3 is the damping coefficient. Their values provide a 
simple and convenient way of characterizing the MOT. 
To find the capture velocity of the MOT, Vc, we con¬ 
sider molecules entering the MOT in the xy-plane and at 
45° to the laser beams, and calculate the fastest speed 
a molecule can have if it is to be captured. We con¬ 
sider MOTs using both the (0-0)A^ni/2 - and (0- 

0)B^E+ - transitions, which we call the ‘A-state 

MOT’ and ‘B-state MOT’ respectively. The lower level is 
the first rotationally excited state of X^E+ (i.e. N = 1), 
while the upper level is the lowest energy state of posi¬ 
tive parity in either the A.^Iii /2 or B^S+ states. Electric 
dipole selection rules for angular momentum and parity 
ensure that these states can only decay back to the = 1 
rotational state of X^E+, meaning that these transitions 
are ‘rotationally closed’. None of the A^n 3/2 - X^E+ 
transitions have this property, which is why we do not 
consider this case. 


III. HYPERFINE STRUCTURE AND ZEEMAN 
SPLITTING 

As we will see, the best scheme for making a MOT 
depends sensitively on the hyperfine structure and Zee- 
man splitting of the ground and excited states. Figure 
shows the relevant energy levels of the X^E+, A^ni /2 and 
B^S+ states, and their Zeeman tuning for magnetic helds 
up to 20 G, calculated using the Zeeman parameters mea¬ 
sured in m- At sufficiently low helds, the Zeeman shifts 
are linear, AE — gp^sMpB, and the values of gp are 
given in Fig.[^ In the X^E+(A^ = 1) state, spin-rotation 
and hyperhne interactions lead to the four components 
shown in Fig. with total angular momenta F =1, 0, 
1, and 2. In the MOT, all four components need to be ad¬ 
dressed. The held-free intervals between these levels are 
calculated using the spin-rotation and hyperhne param¬ 
eters given in HB. The Zeeman tuning is signihcantly 
non-linear for helds above a few Gauss, and there are 
level-crossings above 10 Gauss due to the small energy 
interval between the upper F = 1 and F = 2 levels. 

Figure iRb) shows the much smaller Zeeman tuning 
of the A^i/ 2 ('^ = 1/2) state of positive parity (i.e. 
the e-parity component). The hyperhne splitting of this 
state is unknown. We have attempted to measure it 
by laser spectroscopy using the apparatus described in 
m- The hyperhne splitting of the e-parity and /-parity 
components of the A?IVi/ 2 {J = 1/2) level are, in terms 
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FIG. 1. Zeeman tuning of the states involved in magneto¬ 
optical trapping of CaF. (a) The fonr spin-rotation and hy¬ 
perhne components of the X^S^(w = 0,N — 1) state, (b) 
The two hyperhne components of the A^ni/ 2 (n = 0, J = 1/2) 
state of positive parity (e-parity), (c) The two hyperhne com¬ 
ponents of the B^S^(t; = 0, A = 1) state. The nnmbers in 
brackets are the ^F-factors for each level. 


of the Frosch and Foley parameters, 2/3(/ii/2 — d) and 
2/3 (/ii/2 + d), respectively, where hi /2 = a — {b + c)/2 
and b = bp — c/3. Both hif 2 and d should be small 
and positive (~10MIIz) given the nature of the A^II 
state. Scans over the Pi -h (5i2(l), Qi(0), and ^2(2) 
lines of the (0, 0 )A^ni/ 2 ~Ar^E+ band, which probe the e- 
parity and /-parity components of the A^\iii 2 {J = 1/2) 
level, were recorded at a resolution of 20 MHz FWHM 
and revealed no evidence of excited state hyperfine split¬ 
ting. Similarly, the ■®i?2i(0) and P2(3) lines of the 
(0, 0 )A^n 3/2 — band, which probe the e- and /- 
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parity components of the = 3/2) level, were 

recorded and show no evidence of excited state hyper- 
fine structure. The observed and predicted spectra, and 
associated energy level diagrams, are given in the Supple¬ 
mental Material. These measurements demonstrate that 
/ii/ 2 , /i 3 / 2 (= a-|-(6-|-c)/2) and d are all less than 20 MHz, 
and that the hyperfine interval in the A^Iiii 2 {J = 1/2) 
state is less than 10 MHz. In our model, we set this in¬ 
terval to 4.8 MHz, which is the value determined for the 
hyperfine splitting of the corresponding negative parity 
level [T^ . 

Figure[^c) shows the Zeeman tuning of the = 

0) state. The hyperfine splitting of this state was also 
previously unknown, and so we have attempted to mea¬ 
sure it. Scans over the ^Qi 2 (l) and Pi -|-^ (5ii(l) lines 
of the (0,0)H^I]+ — band were recorded at a res¬ 

olution of 20 MHz FWHM, and they reveal some small 
splittings and shifts due to excited state hyperfine struc¬ 
ture. The observed and predicted spectra, and associated 
energy level diagram, are given in the Supplemental Ma¬ 
terial. From this data we estimate = 20 ± 5 MHz and 
c = 50± 10 MHz. This gives a 20 MHz hyperfine interval 
for the B^E+(iV = 0) state, which is the value we use 
in our model. Because of the small hyperfine splitting 
and rapid Zeeman tuning, the electron and nuclear spin 
decouple at low magnetic fields, and then the states are 
best characterized by the Mj and Mi quantum numbers. 


IV. DUAL-FREQUENCY MOT 

We identify a new mechanism of magneto-optical trap¬ 
ping which avoids the problem of weak trapping associ¬ 
ated with dark states and with a small and which is 
relevant to the rest of our discussion. Consider the simple 
system illustrated in Fig.j^a), which consists of a single 
lower level with FJ = 2, gi = 0.5, an upper level with 
Fu — 1, gu — 0, and two oppositely polarized frequency 
components with detunings (5i and 62 - The wavelength, 
linewidth and mass are set equal to our CaF system. The 
molecule interacts with counter-propagating laser beams 
along each of the coordinate axes, each beam containing 
these two frequency components. The size of the laser 
beams is the same as described above, and the power 
per beam and per component is set to 40 mW. We fix 
61 = —F, and vary ^ 2 . Figure [^b) shows the trapping 
frequency as a function of <52. When 82 is positive, mean¬ 
ing that one component is red-detuned and the other is 
blue-detuned, there is a restoring force. When Mp is pos¬ 
itive the Zeeman shift brings the transition into resonance 
with the red-detuned light, and so we have arranged 
for the red-detuned component of the restoring beam 
to drive AMp = — 1 transitions so that the molecule 
preferentially interacts with that restoring beam (the 
red-detuned component of the anti-restoring beam drives 
AMp = 4-1 transitions). When Mp is negative the tran¬ 
sition comes into resonance with the blue-detuned fre¬ 
quency component, and because it has the opposite po¬ 


larization the molecule again preferentially scatters from 
the restoring beam. The restoring force is maximized 
when <^2 is about O.TF, but remains significant even for 
large positive detunings. There is no trapping when 62 is 
between —0.3F and —F because then the transitions are 
preferentially driven by the anti-restoring beams. How¬ 
ever, for more negative detunings there is once again a 
restoring force. In this case, the restoring beam has one 
component driving AMp = — I transitions, and another 
driving AMp = -|-1 transitions, with the former closer to 
resonance. For the anti-restoring beam, it is the oppo¬ 
site. When Mp is positive, only AMp = —1 transitions 
are driven, and this is far more likely to be driven by 
the restoring beam since it is (relatively) much closer to 
resonance. When Mp is negative, only AMp = 4-1 tran¬ 
sitions are driven, and this is only a little more likely to 
be driven by the anti-restoring beam since the relative 
difference in detuning is not so great. 

Figure [^c) shows how the damping coefficient depends 
on 82 - Cooling occurs when /3 is positive. When 82 is 
negative there is always cooling because both frequency 
components are red-detuned. When <^2 is positive and 
smaller than F there is heating because the blue-detuned 
component is closer to resonance than the red-detuned 
one. When <52 > F there is cooling again, because the 
red-detuned component is the one closer to resonance. 
Together, these plots show that there is both cooling and 
trapping for all values of 82 apart from values between 
0 and —F. It also demonstrates the complexity that can 
arise even in this relatively simple system. 


V. THE A^Ri/z - X2E+ TRANSITION 

We turn now to a model of a CaF MOT using the 
(0,0)A^ni/2 " transition - the A-state MOT. This 

transition has A = 606 nm and F = 27r x 8.3 MHz |12j . 
Four frequency components are used to drive the four 
hyperfine components of the transition, and their po¬ 
larizations in the restoring beam are chosen as shown 
in Fig.|^a), which follows the recommendations of [5] 
for our case where gu is negative. We neglect the small 
branches to higher lying vibrational states. 

Figure [^b) shows how the acceleration of a station¬ 
ary molecule depends on position for various values of 
the power per beam and per frequency component. Here 
we have set the detuning of all components to —F. We 
expect the MOT forces to be weak because of the small 
value of gu^ but we see from the figure that this is not the 
case. This is because the MOT forces are complicated by 
the close spacing of the upper Fi = 1 and Fi = 2 levels 
of the X^I]+ state. This frequency interval is 25 MHz, 
which is only 3F. With our detuning of —F, the field-free 
transition from Fj = 2 is driven (mainly) by two oppo¬ 
sitely polarized components, one detuned by —F and the 
other by 2r. The polarizations of these components are 
exactly as in the model system of Fig.j^a), and Fig.j^b) 
tells us that this configuration produces a trapping force. 
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FIG. 2. (a) Illustration of a model MOT where the lower level has Fi = 2, gi = 0.5, and the upper level has Fu = 1, Qu = O.The 
transition is driven by two oppositely polarized frequency components with detunings of hi and <52, indicated by the dashed 
lines, (b) Trapping frequency versus 52- (c) Damping coefficient versus ^ 2 . 







FIG. 3. GaF MOT operating on the A^ni /2 - transition, (a) The four frequency components of the laser have a 

common detuning, h, from Fi to Tu = 1, as shown by the solid arrows. Dotted arrows show allowed transitions to F^ = 0. 
Values to the left of the arrows give the relative transition strengths for a fixed polarization, summed over and Mi. These 
are independent of the choice of polarization. The polarizations chosen for each component of the restoring beam are indicated 
as CT±. Acceleration versus (b) displacement and (c) speed, for four different values of the power in each MOT beam and each 
frequency component: 5, 15, 40 and 100 mW. The detuning is —F. Acceleration versus (d) displacement and (e) speed, for 
four different values of the detuning: — 0.25r, —0.5r, —F and — 2r. The power is 40 mW. 


Near the centre of the MOT it is this transition which is 
mainly responsible for the restoring force. Further from 
the centre, there is a strong additional trapping force 
coming from the upper F) = 1 transition. At a field of 
about 30 G the Mp = +1 component of this transition 
comes into resonance with the laser component that is 
meant to drive FJ = 2, while the Mp = —1 component 
comes into resonance with the laser component intended 
for F] = 0. Once again, the two laser components are po¬ 
larized correctly so that the transitions from Mp = +1 
and Mp = —1 are both driven by the restoring beam. 
This contribution to the trapping force pushes the turn¬ 
ing point of the acceleration curve further out from the 
centre of the MOT, generating a large MOT capture vol¬ 


ume. We also see from Figure |^b) that the acceleration 
continues to increase as the laser power is increased up 
to 100 mW. 

Figure [^c) shows how the acceleration of a molecule 
at the MOT centre depends on its speed, for the same 
set of laser powers as before. We see that this damping 
force peaks when the speed is near 5 m/s, corresponding 
to a Doppler shift that equals the detuning of —F. There 
is also a secondary peak where the Doppler shift is —dF, 
which is again the result of there being two frequency 
components spaced by SF. This acts to broaden out the 
acceleration curve so that the damping force remains sig¬ 
nificant over a larger range of speeds. 

Figures [^d) and (e) show how the acceleration de- 
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pends on position and speed for various values of the 
detuning, with the power fixed at 40 mW. Bringing the 
detuning close to zero maximizes the trapping force in 
this case, but a negative detuning is required for there to 
be a damping force and so there is a trade-off between 
cooling and trapping. When the detuning is —F and the 
power is 40mW we find = 27r x 160Hz, /? = 5600s“^ 
and Vc = 20 m/s. 

We attribute the strong restoring forces shown in 
Fig-ib ,d) almost entirely to the dual-frequency effect. 
To verify this, we make the intervals between the four 
components of the transition much larger than F, set¬ 
ting the detuning to —F and the power per beam and 
per frequency to 40 mW. Then, we find a very simple 
dependence of the MOT forces on over the range 
0 < gu ^ 0.14, the maximum acceleration is \az,max\ = 
(17000m/s^)|(/„| and the squared trapping frequency is 
ujI = ( 27 rx 282 Hz)^| 5 „|. Note that the sign of Qu dictates 
only the required sign of the field gradient. These rela¬ 
tionships could also be applied to other similar molecules, 
suitably scaled by the ratio of the maximum scattering 
rates. For our present case, = —0.021, and so the 
maximum acceleration is only 357 m/s^ and the trap fre¬ 
quency is 41 Hz. We find similar results if we set the in¬ 
tervals between the four components equal to their true 
values, and also set the polarization handedness of the 
laser components all equal instead of the way they are 
shown in Fig.|^a). These results demonstrate the large 
increase in the MOT forces brought about by the dual¬ 
frequency arrangement. 

It is interesting to consider whether the trapping forces 
can be increased even further by applying the dual fre¬ 
quency method to several of the hyperfine components. 
Figure shows that when = —F and 82 = 2F, the 
trapping and damping forces are both strong, so we can 
aim to arrange this situation for several of the hyper¬ 
fine components. Because the Fi = 2 and upper Fi = 1 
components are spaced by 3F we automatically have this 
ideal situation for the F; = 2 state, as discussed above, 
but this means that we cannot also have it for the up¬ 
per Fi = 1 state. The FJ = 0 state has no Zeeman 
splitting and no dark states, so there is nothing to be 
gained from applying two oppositely polarized frequency 
components. This leaves the lower F; = 1 state. Be¬ 
cause this has a negative 5 -factor, the Zeeman shift of 
the Mp states are opposite to the way they are drawn in 
Fig.ia), and so the polarizations of the two frequency 
components need to be reversed relative to that figure. 
These arguments lead us to add one more frequency to 
address the lower F; = 1 level, detuned by 2F from this 
level and polarized cr—, giving us the set of five frequen¬ 
cies and polarizations illustrated in the inset of Fig.|^ 
Figure |^a) gives the acceleration versus position when 
the laser power is 40 mW in each beam and frequency 
component. The addition of the extra component more 
than doubles the maximum acceleration to 12000 ms“^, 
and the trap frequency increases to = 2tt x 191 Hz. 
Figure |^b) shows the acceleration versus speed for this 
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FIG. 4. Acceleration versus (a) displacement and (b) speed, 
for a CaF MOT operating on the A^ni /2 - transition 

using the set of detunings and polarizations illustrated in the 
inset to (a) and discussed in the text. The power per beam 
and per frequency is 40 mW, and the components are detuned 
by — r apart from the additional component whose detuning 

is 2r. 

case. We find that the additional component has little ef¬ 
fect on the damping coefficient or the maximum damping 
force, but it does narrow down significantly the range of 
velocities where there is damping. In fact, the force has 
the wrong sign for speeds between 8 and 15 m/s because, 
when in the lower F; = 1 level, the molecule is Doppler 
shifted into resonance with the extra frequency compo¬ 
nent of the co-propagating beam. As a result, the capture 
velocity decreases to Vc = 14 m/s. These results suggest 
that a good strategy may be to first load the MOT using 
four frequency components, and then switch on the fifth 
to provide tighter confinement. 

VI. THE B2E+ - X2E+ transition 

Next, we model a MOT using the (0-0)B^E+ - X^E+ 
transition - the B-state MOT. This transition has A = 
531 nm and F = 27r x 6.4 MHz. Population that leaks into 
other vibrational states of can be rapidly pumped 

back into the cooling cycle via the A^ni /2 state. These 
re-pumping transitions contribute very little to the cool¬ 
ing and trapping forces of the MOT, and so we neglect 
them in our model. Four laser frequency components ad¬ 
dress the four transitions from the X^E"*" state, all with 
a common detuning. Their polarizations in the restor¬ 
ing beam are indicated in Fig.j^a), following the rec- 
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(a) 





FIG. 5. CaF MOT operating on the - X^E"*" transition, (a) The four frequency components of the laser have a common 

detuning, <5, from Fi to Fu = 1, as shown by the solid arrows. Dotted arrows show allowed transitions to Fu = 0. Values 
to the left of the arrows give the relative transition strengths for a fixed polarization, summed over and Mi. These are 
independent of the laser polarization. The polarizations chosen for each component of the restoring beam are indicated as (t±. 
(b) Acceleration versus displacement and (c) acceleration versus speed, for four different values of the power in each MOT 
beam and each frequency component: 5, 15, 40 and 100 mW. The detuning is —F. (d) Acceleration versus displacement and 
(e) acceleration versus speed, for four different values of the detuning: —O.SF, —F, —l.SF and —2V. The power is 40 mW. 


ommendations of [8] for a positive We can choose 
whether to excite the transition to = 0 or = 1, 
which are spaced by 20 MHz. The transition strengths 
are given in Fig.j^a). We note that = 0 decays to 
the lower F; = 1 level 95% of the time. If the laser is 
tuned close to resonance with this transition, the combi¬ 
nation of the resonance condition and the stronger tran¬ 
sition strength strongly favors excitation to F^ = 0, and 
so these two levels can form a nearly-closed transition. 
Our model shows that this produces a high scattering 
rate and strong damping, but contributes little to the 
trapping force. Therefore, we choose to drive all F; com¬ 
ponents to F„ = 1. 

Figure [^b) shows the acceleration versus position of 
stationary molecules for four different values of the power 
in each beam and frequency component. The detuning 
is —r. The curves have turning points near the position 
where the largest Zeeman shifts are equal to this detun¬ 
ing. Further out, the acceleration crosses zero because 
the Zeeman shifts become comparable to the upper state 
hyperfine interval and because of level crossings between 
the F; = 1 and 2 levels. The peak acceleration roughly 
triples as the power is increased from 5 to 40 mW, but 
does not increase much for powers exceeding 40 mW. 
This differs from the A-state MOT which requires more 
power to obtain the highest force. For the same set of 
powers, Fig.j^c) shows the acceleration versus speed for 
a molecule at the centre of the MOT. The curves are 
similar to those shown in Fig. I^c) for the A-state MOT, 
except that the range of speeds that give a damping force 


is reduced. The light is red-detuned by 6.4 MHz for exci¬ 
tation to Fu = 1, and so it is blue-detuned by 13.6 MHz 
for excitation to F„ = 0. When the speed is 7 m/s the 
Doppler shift is 13 MHz, bringing the co-propagating light 
into resonance with the transitions to F^ = 0, and thus 
making the net force positive at this speed. This will 
limit the capture velocity of the MOT. 

Figure l^d) shows acceleration versus position for var¬ 
ious detunings, 6 , with the power per beam and fre¬ 
quency component set to 40 mW. As |^| increases the 
turning point of the acceleration curve shifts to larger 
displacements, because this turning point occurs roughly 
where the Zeeman shift matches the detuning. The peak 
restoring force is greatest when the detuning is approxi¬ 
mately —F. A smaller value of |5| reduces this peak force, 
just as it does in a simple one dimensional theory of the 
MOT [T3j. A larger value of |(5| also reduces the peak 
restoring force, which differs from the simple MOT the¬ 
ory, probably because of the different Zeeman tunings 
of the states involved. Figure [^e) shows the accelera¬ 
tion versus speed for the same set of detunings. We note 
that a larger detuning shifts the point where the acceler¬ 
ation crosses zero to higher speeds. When the detuning 
is —F and the power is 40 mW we find = 27r x 254 Hz, 
P = 4600 s“^ and Vc = 10 m/s. Increasing the detuning 
to — 1.5r gives a smaller trap frequency and damping con¬ 
stant, ojz = ‘2tt X 166 Hz, P = 1600 s“^, but the capture 
velocity increases slightly, to 11 m/s, due to the increased 
range of speeds where there is damping. 

When we increase the hyperfine splittings by a 
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FIG. 6. Acceleration versus (a) displacement and (b) speed, 
for a CaF MOT operating on the - X^E"*" transition 

using the set of detunings and polarizations illustrated in the 
inset to (b) and discussed in the text. The power per beam 
and per frequency is 40 mW, and the components are detuned 
by — r apart from the additional component whose detuning 

is 2r. 


factor of 10 from their true values, keeping the laser de¬ 
tunings the same as before, the acceleration curves do 
not change much. This shows that the dual-frequency 
effect discussed in Sec.|IV] does not play a major role in 
this MOT. As for the A-state MOT, it is worth trying 
to engineer a stronger trapping force by using the dual¬ 
frequency effect. Our modelling suggests that a good way 
to do this is to add one extra frequency component with 
tj-l- polarization, detuned by 2r from the upper T) = 1 
level. Figurej^shows the acceleration versus position and 
speed for this arrangement. The detunings and polariza¬ 
tions are illustrated in the inset. With 40 mW in each fre¬ 
quency component, the addition of this extra component 
increases the maximum trapping force by about 50%, and 
the trap frequency increases to = 27r x 278 Hz. The 
damping coefficient and maximum damping force are al¬ 
most unchanged. However, the damping force now has 
the wrong sign for all speeds between 5.7 and 15 m/s, and 
this positive force is larger than before. This is due to the 
Doppler shift bringing the extra frequency component of 
the co-propagating beam into resonance, in addition to 
the the transitions to = 0 being Doppler-shifted into 
resonance as explained above. Consequently, the capture 
velocity decreases to Vc = 8 m/s. 


VII. SUMMARY 

Let us now summarize our findings. We find sim¬ 
ilar trapping forces in the A-state MOT and the B- 
state MOT. For the B-state MOT, when the detuning 
is —F and the power is 40 mW, the maximum acceler¬ 
ation of a stationary molecule is 4800 m/s^, the trap¬ 
ping frequency is uj^ = 27r x 254 Hz, the damping co¬ 
efficient is 13 = 4600 s“^, and the capture velocity is 
Vc = 10 m/s. For the A-state MOT, with the same 
parameters, the maximum acceleration of a stationary 
molecule is 5800m/s^, = 27r x 160Hz, /3 = 5600s“^ 

and Vc = 20 m/s. For the A-state MOT the forces act 
over a wider range of displacements and velocities, which 
is why it has a higher capture velocity. Throughout, we 
have assumed that vibrational re-pump transitions do not 
share the same upper state as the main MOT transition, 
in which case their influence on the MOT forces is small. 
When a re-pump transition does share the same upper- 
state as the MOT transition, we find that the force is 
reduced to about 70% of the values found here. 

We have identified a mechanism of magneto-optical 
trapping that arises when two laser components with op¬ 
posite polarizations and different detunings address the 
same transition. The forces produced are strong even 
when the upper state has little or no Zeeman shift, and 
also when there are dark states amongst the lower lev¬ 
els. In the B-state MOT it is the large value of gu that 
contributes most to the trapping force, with the dual¬ 
frequency effects playing only a minor role. By contrast, 
the trapping force in the A-state MOT comes almost en¬ 
tirely from these dual-frequency effects. These effects 
also play a large role in the SrF MOT that has recently 
been demonstrated igii!. If we turn these effects off 
in our model of the SrF MOT (by artificially increas¬ 
ing the hyperfine intervals) the maximum acceleration 
is reduced by a factor of six, the turning point of the 
acceleration curve is a factor of two closer to the cen¬ 
tre, and the trapping frequency is reduced by a factor 
of two. The dual-frequency effect can be utilized to pro¬ 
duce stronger trapping forces using a carefully designed 
set of laser frequencies. For CaF, adding a single ex¬ 
tra frequency component increases the restoring force by 
100% for the A-state MOT and by 50% for the B-state 
MOT. However, the extra frequency has the effect of re¬ 
ducing the range of velocities where the damping force is 
strong, and this reduces the MOT capture velocity. The 
extra frequency component could be added to tighten the 
confinement once the MOT has already been loaded. 
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